Possibility of spontaneous parity violation in hot QCD 
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We argue that for QCD in the limit of a large number 
of colors, the axial U(l) symmetry of massless quarks is ef- 
fectively restored at the deconfining phase transition. If this 
transition is of second order, metastable states in which par- 
ity is spontaneously broken can appear in the hadronic phase. 
These metastable states have dramatic signatures, including 
enhanced production of r\ and rj mesons, which can decay 
through parity violating decay processes such as r) — > 7r f V°, 
and global parity odd asymmetries for charged pions. 



It may be possible to observe the phase transition(s) 
from hadronic to quark and gluon degrees of freedom 
through the collisions of heavy nuclei at ultrarelativistic 
energies. In the region of central rapidity, the relevant 
phase transitions are those at nonzero temperature; these 
phase transitions can be studied by numerical simulations 
of lattice gauge theory. At present, simulations indicate 
that for three colors coupled to light quarks, there is at 
most one phase transition, controlled by the chiral dy- 
namics of the light quarks The order of the phase 
transition in QCD, in which two flavors are very light, 
and one flavor not too heavy (up, down, and strange), is 
still unsettled. 

The nature of the chiral phase transition depends cru- 
cially upon the dynamics of the axial U(l) symmetry of 
the light quarks ^,|). Notably, for two massless flavors, 
if the axial U(l) symmetry is not restored about the chi- 
ral phase transition, then the transition can be of second 
order; if it is restored, the transition may be driven first 
order by fluctuations. 

There are two approaches to understanding the dy- 
namical breaking of the axial E/(l) symmetry. The first 
assumes that the dominant fluctuations are semiclassical 
instantons §- 0. The second is based upon the large 
N limit of an SU(N) gauge theory ||- [p^| , and assumes 
that the dominant fluctuations are not semiclassical, but 
quantum. 

At zero temperature, both approaches give a reason- 
ably successful phenomenology for the 77' mass and re- 
lated processes. In this Letter we show that these two 
approaches give radically different predictions at nonzero 
temperature. In instanton models of the hadronic vac- 
uum jl, the topological susceptibility is essentially con- 
stant below the phase transition, and only drops off above 



the phase transition. We argue that at large N, the 
topological susceptibility essentially vanishes at the phase 
transition. If the deconfining phase transition is of second 
order, then the axial U(l) symmetry is dynamically re- 
stored as the phase transition is approached from below. 
Under this assumption, using a nonlinear sigma model 
fujl - we show that metastable states with sponta- 
neous parity violation arise in the hadronic phase, and 
would produce striking experimental signatures. 

The large N limit of SU (N) gauge theories is believed 
to be a reasonable approximation even for N — 3 || . We 
assume that confinement holds for all N, with the masses 
of mesons and glueballs of order one as N — > 00; inter- 
actions between mesons and/or glueballs are suppressed 
by powers of 1/N. 

Holding the number of quark flavors fixed as N — > 00, 
at large N the ~ N 2 gluons dominate the ~ N quarks. 
Taking the degeneracy of hadronic bound states to be of 
order one, the gluonic free energy changes from ~ N° 
in the hadronic phase, to ~ N 2 in the deconfined phase. 
Thus the gluonic part of the free energy can be used to 
define the temperature of the transition, at T = T<j ~ N° 
assume that any other transitions in 
the theory also occur at T<j. Given the huge change in 
the free energy, any other possibility seems baroque, at 
best. 

In the pure glue theory, the topological susceptibility 
Xym{T) = d 2 F{0,T)/d0 2 = J d 4 x (Q(x)Q(0)), where 
F(6,T) is the free energy, and the 9 parameter is con- 
jugate to the integral of the topological charge density, 
Q{x) = (g 2 /32Tr 2 )tr(G al3 G a P). At zero temperature, the 
free energy reduces to the energy, F(6, 0) = E{6). 

Because Q(x) — d a K a , where K a is the (gauge vari- 
ant) topological current, Xym{T) vanishes order by or- 
der in perturbation theory. At zero temperature, Witten 
suggested that quantum fluctuations generate a nonzero 
value, \ym(0) ~ A° |]]. At high temperature, the theory 
is asymptotically free and so weakly coupled, with elec- 
tric fluctuations suppressed by Debye screening. Thus 
at high temperature, \ym{T) is unequivocally calcula- 
ble by semiclassical means, using instantons ^ffl. With 
g 2 the gauge coupling, the instanton action is 8ir 2 /g 2 ; as 
g 2 N is held fixed when N — > 00, \ym{T) ~ exp(-aN), 
a = 8ir 2 /(g 2 N). This naive picture was verified by Af- 
fleck in a soluble asymptotically free theory, the CP N 
model in 1 + 1 dimensions ||. 

Thus for gauge theories, Xym(T) changes from ~ N° 
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at low temperatures, to ~ exp(-aN) — which at large 
N is essentially zero — at high temperature. Appealing 
to simplicity, we assume that this change happens at the 
deconfining transition, X YM (T) ~ for T > T d §,[[§. 

How \ym{T) vanishes as T — > T7 depends upon the 
order of the deconfining phase transition. If the decon- 
fining phase transition is of first order, then since all in- 
teractions at large N are suppressed by 1/N, presumably 
Xym{T) = Aym(O) for T< T d . 

While there is some evidence that the deconfining 
phase transition is of first order for all TV > 4, this con- 
clusion may be premature Q . Following the conjecture 
of Jl9| , we henceforth assume that the deconfining phase 
transition is of second order at large N . This requires 
that the phase transition is driven by a Hagedorn spec- 
trum |^7| . Adding Nj > 2 flavors of massless quarks, 
we assume that chiral symmetry is spontaneously broken 
at zero temperature ||, and that the deconfining phase 
transition forces chiral symmetry restoration at T = T d . 
Since we shall argue that A(T) — > as T — > T^~, the 
chiral transition is driven first order by fluctuations . 
At large N, however, since the mesonic couplings start 
out small, ~ 1/N, fluctuations can only drive the chi- 
ral transition first order within a narrow critical region, 
~ 1/N. Thus the latent heat for the chiral transition is 
~ N°, and relative to the ~ N 2 for the gluons in the de- 
confining phase, it is a very weakly first order transition. 
Consequently, we let quantities associated with the chi- 
ral transition vary with temperature, but only as in mean 
field theory. Notably, the pion decay constant f„, which 
is ~ VN, is assumed to decrease as f n (T) ~ (T d — T) 1 / 2 
when T — > T d . 

While the pure glue theory depends upon 9, the ad- 
dition of massless quarks must cancel any 9 dependence. 
Witten showed that this happens at large N by the ap- 
pearance of a light meson, the rf ||. We generalize this 
to nonzero temperature, to estimate how Xym(T) and 
the n' mass vanish as T — > TT . 

At large N, at zero temperature any gauge invariant 
correlator is saturated by the exchange of single glueballs 
and mesons || . Normally this changes in a thermal bath, 
due to scattering off states in the thermal distribution. 
The hadronic phase of large N QCD, however, is "cold": 
chiral symmetry is restored not at a scale set by the pion 
decay constant, /„• ~ \/~N (as in, say, the sigma model at 
large N) , but at a much lower temperature, given by the 
deconfining transition, T d ~ N°. Thus we can use the 
same type of arguments as at zero temperature, simply 
allowing any quantity which enters to be temperature 
dependent. 

We start by following Veneziano , and define A, ; > as 
the form factor between the topological current and the 
rf meson, (0\K a \n') = i(y/Nf /N)p a X n > (T), with p a the 
momentum of the rj' meson. This form factor is precisely 
analogous to the coupling of ir° to two photons. Following 
, to one loop order in a constituent quark model, the 



(anomalous) coupling of the rj' to two, or indeed any 
finite number of gluons, vanishes as chiral symmetry is 
restored, like A V (T) - f„(T) ~ (T d - T) 1 ' 2 as T -> T^ . 
Using Veneziano's relation, m 2 ,(T) — Nf X 2 , (T) / 'N , we 
find the rf mass vanishes as m 2 ,(T) ~ (T d — T) when 

We now use Witten's formula B for the rf mass, 
m 2 ,(T) = AN f X Y M{T)/f 2 (T). This relation shows that 
Xym{T) and Jtt(T) must vanish at the same point, but 
not much else. Since Veneziano's formula tells us how 
m 2 , (T) vanishes, though, we can use this to deduce how 
the free energy depends upon 9 about T d : 

F(9,T)~(l + c9 2 )(T d -T) 2 - a , T^T rf -, (1) 

for some positive constant c, \9\ < it. That the 9 de- 
pendence is only quadratic is characteristic of large N 
§. Then X YM (T) = d 2 F{9,T)/d9 2 - (T d - T) 2 ~ a and 
m 2 ,(T) - (T d - Tf- a as T -> T7 . Because the critical 
exponent a^O, this does not quite agree with our esti- 
mate using Veneziano's formula. We trust (|j), since the 
calculations of |2(| are only at one loop, and so basically 
mean field. Even so, as a ~ —.013 O, this difference 
is small, and does not alter the conclusion that Xym(T) 
and m^^T) vanish at the phase transition, with T d and 
a independent of 9. 

Previously, Affleck p| and also Davis and Mathieson 
H] argued that Xym(T) vanishes when T > T d \ our con- 
tribution is to estimate how it vanishes as T — > if the 
deconfining phase transition is of second order. If the de- 
confining transition is of first order, then as the hadronic 
phase is cold, presumably Xym{T) = A-km(Q) for T < T d , 
at which point it drops discontinuously to zero. (For an 
alternate view, with T d ^T x , see Meggiolaro JOf.) 

An effective nonlinear sigma model which incorporates 
the breaking of the axial U(l) symmetry can be con- 
structed jijj- U|. With U a U(N f ) matrix satisfying 
TftU = 1, the potential for U is 

f 2 

V(U) = ^ (tr (M(U + C/t)) _ a (tr In U - 9) 2 ) ; (2) 

M is the quark mass matrix. When M = 0, m 2 , ~ a, so 
a - X 2 ,/N. (Our a = a/N in 0- @.) 

Instanton processes are often modeled using a linear 
sigma model with a field $, by introducing a term ~ 
e ie det(§) 0; this term is well behaved in both the 
low and high temperature phases. In contrast, in (||) the 
term ~ (tr In U) 2 , which breaks the axial U(l) symmetry, 
only makes sense in the low temperature phase, and is 
singular if the vacuum expectation value (v.e.v.) of U 
vanishes. At large N, however, everything fits together: 
since a(T) — > as T — > TJ , there is simply no such 
term in the high temperature phase. Moreover, at large 
A, even at zero temperature ~ det(<&) must be dropped, 
since it is inconsistent with the 9 dependence p3]. 
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Taking My = (J^S^ , any v.e.v of U can be assumed to 
be diagonal, Uij = e^S 1 ^; then 

V(4>i) = coaM + l^h-oA , (3) 

which is minimized for sin(<pi)+a(^2 4>i — 9) = . Note 
that as c/>i arises from tr In U, it is defined modulo 2ir. 

All of the parameters in (|^) are temperature depen- 
dent. When a = 0, the Goldstone boson masses fi 2 ~ 
irii(qq)/ f 2 , with rrii the current quark masses, and (qq) 
the chiral order parameter. When M ^ 0, there is no 
true critical point, but we can use mean field theory to 
estimate that as T — ► TT , / w and a decrease, while the 
Hi all uniformly increase: / 2 (T) - a{T) ~ (T d - T), 
(qq) ~ (T d - T) 1 / 2 , and M |(T) ~ m,/^ - T) 1 / 3 . The 
solutions are independent of /„■, and depend only upon 
the ratio a//xf . In mean field theory this ratio is indepen- 
dent of flavor, and scales as a{T)/u 2 (T) ~ (T d - T) 3 / 2 . 

Several authors have studied how the v.e.v's for the 
(j>i change at zero temperature when 9^0 
Instead, we consider 9 = 0, and follow Witten ]13| to 
investigate metastable solutions at small a. Related 
metastable states have been discussed by Shifman p2]] . 

For a single flavor, the vacua are at <p = 0, ±27r, etc. By 
balancing [i 2 sin(<p) against a0, however, for small a//i 2 
it is easy to show that there are other solutions with 
<p =^ 0. These solutions have higher potential energy, and 
so are local but not global minima. Numerically, we find 
that the first metastable state occurs when a < a cr , with 
o-cr/ IJ 2 ~ .217, and (p c ~ 4.493; the (p field is massless 
about <p c . As a — > 0, </> — > 2ix, which is equivalent to 
cp = 0. There is an infinite tower of metastable states; we 
only consider that with lowest energy. 

These metastable states are like regions with nonzero 9, 
and so spontaneously break CP symmetry. Under charge 
conjugation, (p — > +<p, while under parity, (p — > —<p. Al- 
though there is a solution at —<f>, when —<f> does not differ 
from +<p by a shift of 2n, parity is spontaneously vio- 
lated. This does not conflict with Vafa and Witten [El], 
who showed that at 9 = 0, parity is not spontaneously 
violated in the QCD vacuum. Their theorem generalizes 
to the thermodynamic minimum at nonzero temperature 
(although not at nonzero quark density), but it does not 
constrain metastable states. 

The appearance of metastable states when Nt > 2 
is somewhat subtle. To illustrate the basic point, we 
consider a two flavor model in which u 2 = and /i 2 , = /i 2 . 
If the two flavors decoupled |U5|, one might guess the 
solution (pi — and (pi as for one flavor. The equation 
of motion for (pi, however, forces (pi + (pi — 0: when 
fi\ = 0, for any value of a there is only the trivial solution, 
4>i = 02 = 0, modulo 2n. 

This example shows that there are no metastable states 
if any quark mass vanishes. This is natural: after all, 
there is also no 9 dependence when any quark mass 



vanishes, and these metastable states are similar to re- 
gions with nonzero 9. Analogously, when all quarks have 
nonzero mass, metastable states only appear when a is 
small relative to the lightest quark mass. Thus in QCD, 
for metastable states to occur a must be small relative 
not to the strange quark mass, but to the up and down 
quark masses. For this reason, a(T) must become very 
small near the phase transition. 

The potential of (|) can be used to obtain a qualitative 
estimate. Let m u , mj, and m s be the masses of the 
up, down, and strange quarks. The charged pions and 
the kaons are unaffected by the anomaly; with ti 2 ± ~ 
m u + md and m 2 K ~ m Ui d + m s , and assuming that m u = 
m d /2, m n ± = 140 MeV and m K ~ 496 MeV give /i 2 = 
(lUAfe^) 2 , i4 = (161 MeV) 2 , and //§ = (687 MeV) 2 . 
We take Veneziano's value of a = (492 MeV) 2 , and 
numerically diagonalize the mass matrix in (JsJ) to obtain 
TO^o - 139 MeV, m v - 501 MeV, and m, r - 983 MeV. 
This is reasonably close to the experimental values of 
m v - 548 MeV, and m v > ~ 958 MeV. 

Taking the ratios of the fi 2 as fixed, and varying a//i 2 , 
we studied numerically the appearance of the lowest en- 
ergy metastable state. For the sake of discussion, we take 
the zero temperature /i 2 ; thus only the ratios of masses 
are believable, with all true masses larger by some uni- 
form factor, ~ /j 2 (T)//j 2 (0) - (T d - T)~ x l 2 . The masses 
of the 7T°, 77, and rj' are read off by diagonalizing the 
mass squared matrix obtained from (|3|). The masses of 
the charged pions are m 2 ± = a\cos((pi) + /jjcos^)- We 
ignore changes in the kaon masses; as (p3 is small, their 
masses do not change much. 

We find that there is a metastable solution when 
a/i^ 2 < .2467, but it is unstable in the ir° direction un- 
less a < a cr , a cr / ' n\ ~ .2403. At a cr , <pi ~ 4.47, <pi ~ 
— .524, and 03 ~ —.028; the 7r° is massless at a cr , while 
Tn 7r ± ~ 106 MeV, m n - 150 MeV, and m v > - 687 MeV. 
As a — > 0, the metastable state becomes equivalent to 
the vacuum, as <pi — > 2ir, <pi and (p3 — > 0. At a = 0, 
TO^o ~ HAMeV, m^± ~ 140 MeV, m, - 161 MeV, and 
rriTji ~ 687 MeV. In the thin wall approximation |23], the 
decay rate of the metastable state is T « exp(— F c /T), 
where F c ~ (32 % /2/3)(/ J ?/ 2 / a 2 ). 

Putting in the zero temperature values, in order for 
metastable states to occur, near T d the ratio of a/ fi 2 must 
be about 1% of its value at zero temperature. It is not 
clear if this is possible in QCD, but of course this esti- 
mate is manifestly model dependent. Since F c ~ 1/a 2 , 
at small a the metastable states live a very long time; 
thus in heavy ion collisions, metastable states do not de- 
cay by bubble nucleation; instead, as the hot phase cools, 
the value of a changes dynamically, and the metastable 
state rolls smoothly into the true vacuum. 

When a becomes very small, there are several features 
common to both the ground state and the metastable 
states. First, the neutral Goldstone bosons are eigen- 
states not of 577(3), but of flavor ||,|l]: at a = 0, 
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7T° ~ uit, 77 ~ dd, and 77' ~ ss. This generates maxi- 
mal isospin violation: the neutral pion is lighter than the 
charged pions, and so produced more readily. This effect 
is much stronger for the metastable states, since the tt°'s 
are massless at a cr , and so very light for a ~ a cr . Simi- 
larly, the 77 and 77' also become light in both phases; this is 
especially true for the 77, as it sheds all of its strangeness. 
It is not clear how much lighter the 77' becomes, given 
the overall increasing mass scale of /if (T). Light rj and 
77' mesons are produced more readily ]24| , and can be 
observed either directly, through 77 decays p4j] , or indi- 
rectly, through pion Bose- Einstein correlations |^5j . 

There are two types of experimental signatures special 
to the formation of a parity violating phase. The first 
is that decays normally forbidden by parity are allowed 
(for a related phenomenon, see J26|). Kinematically, rj — > 
7r + 7r~ is not allowed, but 77 — > 7t°7t° is. The processes 
77' — > 7r + 7r~ and 77' — > 7r°7r° are also allowed; however, as 
the 77' is almost pure ss, this is suppressed by ~ m u /m s . 

There are also global variables which are sensitive to 
the dynamics of a parity violating phase. It can be shown 
that the interactions of charged pions differ if there are 
regions with <f>\ and <f>2 7^ 0, which change in either space 
or time. This is similar to the propagation of charged 
particles in a background magnetic field: an e + e~ pair, 
produced back to back, are both deflected in the same 
direction by a magnetic field. A parity odd asymmetry 
could be observed by summing over all tt + tt~ pairs in a 
given event, 
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z is the beam axis of the collision, and P^± are the pion 
momenta. V is like handedness in jet physics j27j. 

These metastable domains might be of cosmological 
interest. A region with fa 7^ implies a s G a pG al3 ~ 
^ym(T) 4>i [0; likewise, the coupling to electromag- 
netism should also generate aE-B ~ a s G a pG al3 . Thus if 
the entire universe fell into such a metastable domain, it 
would generate a nonzero value for a cosmological mag- 
netic field at the time of the QCD phase transition p8[ . 

To summarize, in the limit of large N the topologi- 
cal susceptibility (essentially) vanishes in the deconfincd 
phase. If the deconfining transition is of first order, then 
the susceptibility is constant in the hadronic phase; if of 
second order, the susceptibility vanishes as T — > Tj, (1). 
If the latter happens, metastable states in which parity 
is spontaneously broken can appear, although in a non- 
linear sigma model, one must be very close to the phase 
transition for them to occur. It is not clear if the large 
N expansion is a good guide to this physics when N = 3; 
certainly at N = 3, the susceptibility will be nonzero 
in the high temperature phase. Nevertheless, the large 
N expansion provides a qualitative guide against which 
other models can be tested. 
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